A Critical Mach Number for Electron Injection in Collisionless Shocks 
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Electron acceleration in collisionless shocks with arbitrary magnetic field orientations is discussed. 
It is shown that the injection of thermal electrons into diffusive shock acceleration process is achieved 
by an electron beam with a loss-cone in velocity space that is reflected back upstream from the 
shock through shock drift acceleration mechanism. The electron beam is able to excite whistler 
waves which can scatter the energetic electrons themselves when the Alfven Mach number of the 
shock is sufficiently high. A critical Mach number for the electron injection is obtained as a function 
of upstream parameters. The application to supernova remnant shocks is discussed. 

PACS numbers: 95.30.Qd, 96.50.Fm, 96.50. Pw, 98.38.Mz, 98.70.Sa 



It is widely believed that a large part of cosmic rays 
(CRs) in our galaxy are produced by supernova remnant 
(SNR) shocks. The diffusive shock acceleration (DSA) 
model, which is based on the idea of particle scattering 
by plasma waves, has been considered as the most plausi- 
ble mechanism of CR acceleration at collisionless shocks 
Since ions can easily be scattered either by ambient 
or self-generated magnetohydrodynamic (MHD) waves, 
the ion acceleration (or injection) is relatively easier to 
understand @). On the other hand, resonant scatterings 
of thermal electrons are usually less efficient because of 
their small gyroradii, and thus, mildly relativistic ener- 
gies are needed for the resonance with MHD turbulence 
in the typical interstellar medium. This, the so-called 
injection problem, has been a long-standing unresolved 
issue in the DSA theory. Indeed, it is known from in situ 
observations in the heliosphere that the gradual increase 
of energetic particles toward the shock, as predicted by 
the DSA theory, is the common feature for ions, while 
it is not for electrons. The accelerated electrons (typi- 
cally with energies < lOkeV) are usually confined in a 
much smaller region close to the shock 0, 0] except for 
rare occasions [1, 0]. These observations in the helio- 
sphere support the theoretical expectation that the in- 
jection of electrons is far less efficient than that of ions. 
On the other hand, there is strong evidence for the pres- 
ence of ultra-relativistic electrons in young SNR shocks 
. The most striking difference between the heliosphere 
and other astrophysical environments seems to be the 
difference in Alfven Mach numbers Ma = V s /va, where 
V s , and va are the shock, and the Alfven velocity, re- 
spectively. The solar wind at 1 AU has an average Mach 
number of ~ 5—10, while it can be much higher at SNR 
shocks. The purpose of this paper is to show that there 
exists a critical Mach number above which the electron 
injection can naturally be explained. 

The injection of electrons requires either a mechanism 
by which thermal electrons are accelerated to mildly rel- 
ativistic energies, or the generation of high-frequency 
whistler waves propagating toward the shock which can 



scatter sub-relativistic electrons. Levinson [8j proposed 
the injection process through the excitation of high- 
frequency (oblique) whistler waves by the anisotropy of 
pre-existing CR electrons. A different approach is to 
consider the pre-acceleration of thermal electrons to the 
injection threshold through plasma microinstabilities in 
the thin shock layer. One possible scenario is to con- 
sider lower-hybrid waves excited by the reflected ions in 
quasi-perpendicular shocks (0B n ^ 45, where the shock 
angle Obu is the angle between the shock normal and 
the upstream magnetic field), while it requires a rather 
high shock velocity (a few per cent of the speed of light), 
as well as large magnetic field for the injection An- 
other mechanism working in weakly magnetized plasmas 
is the shock surfing acceleration (SSA) [l(J which 
was recently shown to be efficient even in two dimen- 
[la ). Amano and Hoshino [3] have recently 
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demonstrated using one-dimensional particle-in-cell sim- 
ulations that suprathermal electrons accelerated by the 
SSA are subject to further energization by the shock drift 
acceleration (SDA). The process can actually be a possi- 
ble solution to the electron injection problem because the 
accelerated electrons can self-generate Alfven waves. We 
begin our discussion by deriving the condition required 
for the self-generation of Alfven waves. 

The SDA can be understood as an adiabatic mirror 
reflection process in the de Hoffman- Teller frame (HTF) 
where the motional electric field vanishes. The beam 
velocity of the reflected population measured in the up- 
stream rest frame is given by 14 ~ 2V S / cos 6>_b„. We 
consider electromagnetic instabilities excited by the re- 
flected electron beam. In the following, the mass, charge, 
density and temperature of particle species s are respec- 
tively given by m s , q s , n s , T s . We then define the plasma 
frequency, gyrofrequency, thermal velocity, Alfven veloc- 
ity, and plasma beta as follows: oj ps — y/ '47m s q 2 /m s , 
O cs = q s B/m s c, v s = yj2T s /m s , va, s = B / yJAiin s m s , 
Ps = v s/ v a si where B is the ambient magnetic field 
strength. Note the definition of the gyrofrequency in- 
cludes the sign of charge (i.e., £l ce < 0). Figure Q] 
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FIG. 1: Schematic dispersion diagram for circularly polar- 
ized electromagnetic waves in an electron-ion plasma. Posi- 
tive (negative) frequency corresponds to the right-hand (left- 
hand) polarization. The cyclotron resonance condition ui = 
kVb — fl C e(Vb < 0) is also shown. Waves in the shaded regions 
are strongly damped by the cyclotron damping of thermal 
plasma. 



schematically shows the dispersion relation of circularly 
polarized electromagnetic waves propagating parallel to 
the magnetic field in an electron-ion plasma. Note that 
the positive k direction is pointing toward the shock (thus 
Vi, < 0). The beam strongly interacts with the back- 
ground plasma only at the intersection points between 
the resonance condition lu = kVb — O ce and the nor- 
mal mode of the plasma (A and B in the figure). If 
one considers a cold beam, i.e., 8 function in velocity 
space f{v\\ , v±) oc S(vh — Vb)S(v±) , v± are the velocity 
components parallel, and perpendicular to the ambient 
magnetic field, respectively), the excitation of whistler 
waves at the point A is prohibited because otherwise the 
momentum conservation law is violated. On the other 
hand, the wave growth at the point B should overcome 
the cyclotron damping of thermal ions, which is strong 
at short wavelength kc/u p i > 1 (where c is the speed of 
light). Therefore, the excitation of Alfven waves requires 
that the reflected electron beam interacts with the nor- 
mal mode at kc/ui p i 1. Assuming \lo\ <C \fl r 
obtain the condition 11^ 



M A > 



cos 9 Bn mi 



(1) 



The required Mach number is fairly high especially at 
quasi-parallel shocks. Although some young SNRs may 
satisfy the condition, it seems to be too stringent to ex- 
plain available astrophysical observations. 

We now want to relax the requirement by exploring 
the possibility of the whistler wave excitation at the 
point A. It is actually possible by considering an elec- 
tron beam with a loss cone in velocity space. The loss- 
cone distribution is a natural consequence of the mir- 



ror reflection because particles having large pitch angles 
are preferentially reflected by the magnetic mirror. It 
can be unstable against the excitation of whistler waves, 
however, the instability is tend to be suppressed by the 
cyclotron damping of thermal electrons. The damp- 
ing of whistler waves is significant for waves satisfying 
—kv e — fl ce < to < +kv e — fl ce (the upper shaded region 
in Figured]). Therefore, the wave generation is expected 
only when the beam velocity is faster than the electron 
thermal velocity \Vb\ > v e . The condition can be rewrit- 
ten as 



M A > C ^ x ^ = Mf- 
2 \ m e 



(2) 



Since the required Mach number is now proportional to 
y/ rrii/m e , it is smaller than equation ([!]) by a factor of 
~ 43/v^e, thereby greatly relaxes the requirement for 
plasmas with B e ~ 1. 

The above analytical expression gives only a rough es- 
timate. While it is useful, we want to determine the con- 
dition more precisely. To do so, we need to model a veloc- 
ity distribution function of the reflected electron beam. 
We here adopt a cold ring-beam distribution of the form 
f(v\\,vj_) = n b /(2wV r )S(v\\ - Vb)5(vj_ - V r ) where n b , V r 
are the beam density, and "ring velocity", respectively. 
Although it may not necessarily be a realistic distribu- 
tion, the essential physics is retained for our purpose. 
The cold plasma dispersion relation for right-hand circu- 
larly polarized electromagnetic waves in the presence of 
the cold ring-beam electrons can be written as 
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where 77 = rih/n, <C 1 is the density of the beam nor- 
malized to the total density. We seek an approximated 
unstable solution of the dispersion relation in the low 
frequency regime (kc/uj 3> 1). It is easy to understand 
that the contribution from the beam becomes significant 
when uj ~ loq = kVb — C e- Thus, we retain only the 
term proportional to l/(w — uiq) 2 for the beam. We de- 
fine fco as the wavenumber at U) = ujq on the whistler 
mode branch in the absence of the beam. We then write 
lu = luq + 5u> and expand the dispersion relation in pow- 
ers of Slo under the assumption \5ui\ -C \(jJo\- The growth 
rate 7 gr0 w = Im(Sui) is then obtained as 
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where x = -(u + fl ci ) 2 (io + n ce ) 2 / (n cl (fl ci (io Q + ri ce ) 2 + 
O ce (wo + n r i) 2 }mi I m e ) i s a positive numerical factor of 
order unity [see, e.g.,[f5|. It is important to note that 
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the growth rate is proportional to V r , in addition to 
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the well-known dependence for beam-plasma instabilities 
77 1 / 3 . The wave growth, in reality, competes with the the 
cyclotron damping by a thermal plasma. The cyclotron 
damping rate 7damp for a wave (wo,fco) may be written 
as 

7damp <? s. - n s m, fl ci - ( 1+ ) 

fl ci uodD r /du)\ u = Uo v\ m m s k v s 

where D r is the real part of D. For calculation of the 
damping rate, the contribution from the beam compo- 
nent is ignored. We assume that the ion temperature is 
the same as the electrons — j3 e ), and consider the 
damping rate as a function of j3 e for a given {u)Q,ka). 
Since (wo,fco) depends only on V&, the total growth rate 
7tot = 7grow + 7dam P is solely determined by V b , V r , /3 e , 
and rj. In the present paper, the beam density rj is con- 
sidered as a free parameter, on which the dependence of 
the result should be investigated. 

For our purpose, it is sufficient to write the ring veloc- 
ity V r by using upstream quantities. It may be estimated 
as the minimum perpendicular velocity required for the 
mirror reflection 

v 2 ^ (v b 2 + ^4>vAtzn 2 e lc , 

where cf> — 2e0 HTF /miV s 2 , and 9i c are the normalized 
cross-shock electrostatic potential measured in the HTF, 
and the loss-cone angle, respectively. Since the potential 
is estimated to be 0~ 0.1—0.2 from in situ observations 
[ltij . a fixed value of = 0.2 is used. The loss-cone an- 
gle 9i c is defined as tan 2 9i c = B /(B max — B ) using the 
magnetic field strength in the upstream B and at the 
maximum B max . It is worth noting that the compressed 
transverse magnetic field reaches its maximum at the so- 
called overshoot region, where the compression ratio r 
is greater than expected from Rankine-Hugoniot rela- 
tions. Since the foreshock region of quasi-parallel shocks 
are always accompanied by large amplitude MHD tur- 
bulence, it is natural to assume the presence of fluctu- 
ating transverse magnetic fields 5B ~ Bo ahead of the 
shock. Hence, we may estimate the maximum trans- 
verse magnetic field as ~ SB + rsm9 Bn B . By assum- 
ing SB = Bq, we obtain the maximum magnetic field 
strength {B max /B ) 2 = cos 2 9 Bn + (r sin 9 Bn + 1) 2 . Thus, 
we obtain a finite compression ratio even at a purely par- 
allel shock (9 Bn — 0). We assume r — 6 in the following 
discussion. 

Since both V& and V r are written in terms of Ma and 
9 Brll the total growth rate 7tot is now a function of Ma, 
9 Bn , and /3 e for a given r\. We define a critical Mach num- 
ber M™ J * for the electron injection as a numerical solu- 
tion of 7tot = 0. Figure [2] shows the solutions as functions 
of 9 Bn for 77 = 10~ 3 , 10~ 4 , 10~ 5 , where a fixed f3 e = 1 is 
used. For comparison, a solution with a constant tan 9i c 
assuming B max — tBq is also shown by a dashed line 
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FIG. 2: Critical Mach number for the electron injection as a 
function 0Bn for the case f3 e = 1. 77 = ni/rii is the normal- 
ized beam density. The dashed line corresponds to the case 
with a constant #; c . The black solid line shows the analytical 
expression given in equation (J2J) for reference. 

for the case of r; = 10 -4 . These two curves are almost 
identical in the quasi-perpendicular regime, while the de- 
viation becomes larger at quasi-parallel shocks. We find 
that the result only weakly depends on the choice of 0, 9i c 
at 9 Bn > 45 for the parameter range of our interest. On 
the other hand, the result may depend on the assumption 
at a quasi-parallel shock due to its intrinsic complexity. 
For instance, our model expects large ratios of V r /Vb for 
quasi-parallel shocks (e.g., K/H — 4.7 at 9 Bn = 30) 
compared to quasi-perpendicular shocks (V r /Vb < 1), 
which makes it easier to excite whistler waves. Although 
this trend itself is naturally expected, specific quantities 
will certainly depend on the shock structure. Therefore, 
we think that our estimate in the quasi-parallel regime 
includes an uncertainty factor of order unity, which can 
be studied only by using self-consistent numerical simu- 
lations. 

The above result well explains a recent statistical anal- 
ysis of the Earth's bow shock crossings observed by the 
Geotail satellite Q. Oka et al. 4] claimed that the 
power-law index of accelerated electrons T clearly de- 
pends on the so-called whistler critical Mach number 
M whisticr = ( a /2) C os 9 Bn ^/mi/m e , which is defined as 
the critical point above which whistler waves cannot 
propagate upstream. Note that a is a numerical fac- 
tor of order unity, and they used a = \/27/16 ~ 1.3. 
It is easy to recognize that M][ at j g c \ ose ^ j n 
plasmas with /3 e ~ 1 , which is statistically valid in the av- 
eraged solar wind. We see from their Figure 4(b) that T 
flattens when Mach numbers go slightly above MJ hlstler . 
In a parameter regime relevant to their observations (i.e., 
10~ 5 < ^ < 10~ 3 , P e ~ 1, and 9 Bn > 60), we find roughly 
M^*/M^ - 2 or slightly below. The result is consistent 
with the observation in that the critical Mach number for 
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the efficient electron acceleration is larger than j\^whisticr 
by a factor of ~ 2. Note that our extensive parameter 
survey shows that roughly M_^ J */M™ J ^1—3 for a more 
wide range of parameters. 

Since the growth rate of whistler waves well above the 
threshold is much larger than the ion gyrofrequency (e.g., 
Ttot/Aci 21 for Ma = 10, /3 e = 1, 6 Bn = 80, and 
rj = 10~ 4 ), the spatial scale of the wave growth is much 
smaller than the convective gyroradius of ions Vg/jtot *C 
V s /Q c i, which represents the typical width of the quasi- 
perpendicular shock transition region. Therefore, we ex- 
pect that the electron acceleration occurs within the thin 
shock transition region, which again well explains the ob- 
served characteristics of energetic electrons. We should 
note that loss-cone type distributions are frequently ob- 
served in the foreshock region and thought to be a result 
of the mirror reflection process [I3|. Therefore, the as- 
sumption of loss-cone distributions seems to be reason- 
able at first glance. On the other hand, the presence of 
a loss cone in the upstream region means the absence of 
whistler-wave generation and resulting pitch-angle scat- 
tering, which is most likely due to lower Mach num- 
bers. Careful reanalysis of the electron measurements 
as a function of Mach numbers within or immediate up- 
stream of the shock is needed to further ensure the appli- 
cability of the present theory. It is important to mention 
that Gosling et al. Q reported that suprathermal elec- 
trons are observed primarily around quasi-perpendicular 
shocks, with an exception of an unusually high Mach 
number solar wind Ma ~ 41. According to the authors, 
the shock (#b„ ~ 45) was quasi-parallel in character but 
has an enhanced suprathermal electron flux, which is also 
consistent with our prediction. 

So far, we have confirmed that the present theory 
is strongly supported by in situ measurements of the 
bow shock. It is interesting that the critical parame- 
ter Ma/VPc is independent of both the ambient density 
and the magnetic field strength. Our result indicates that 
the required condition for the electron injection may be 
written by using typical SNR shock parameters: 



M A 



68 



3000 km/sy V 10eV 



/ T 



-1/2 



> 30 



(3) 



regardless of the upstream magnetic field directions 
(#B n ). Therefore, we conclude that the injection and 
subsequent acceleration of electrons through DSA pro- 
cess will operate in young SNR shocks. This gives, for 
the first time, the most natural explanation for the dif- 
ference in the electron acceleration efficiencies observed 
in different environments. 

In the present study, we have performed only the linear 
instability analysis. However, it is known that whistler 
waves undergo both forward and inverse cascade pro- 
cesses (llj . which may enhance the efficiency of injec- 
tion and further acceleration to relativistic energies. We 



should note that a quantitative estimate of the injection 
rate r/ ~ 10 -4 is given by [l4[ with some simplifying as- 
sumptions. Although it is possible to consider nonlinear 
processes affect the estimate quantitatively, the injection 
process itself will not be modified as it has only a weak 
dependence on the injection rate. Therefore, we think 
that the present mechanism for the electron injection is 
robust, and there is no reason why the process should 
not work in astrophysical shocks. 

One of the most important remaining problems, when 
it applies to SNRs, is probably the effects of CR back 
reaction. It is believed that CRs amplify the upstream 
magnetic field up to ~ 0.1—1 mG as suggested by X-ray 
observations of some young SNRs. While the condition 
given by equation <j3j> is independent of the magnetic field 
strength, other effects such as pre-shock electron heat- 
ing and deceleration of the upstream plasma, if strong, 
might become important. Although we think the injec- 
tion mechanism itself will remain unchanged because of 
the locality of the process, it is interesting to investigate 
the consequence of such nonlinear shock behaviors to the 
electron injection efficiency. 
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